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Abstract A search for a fermiophobic Higgs boson us- 
ing diphoton events produced in proton-proton colli- 
sions at a centre-of-mass energy of -/s = 7 TeV is per- 
formed using data corresponding to an integrated lumi- 
nosity of 4.9 fb _1 collected by the ATLAS experiment 
at the Large Hadron Collider. A specific benchmark 
model is considered where all the fermion couplings to 
the Higgs boson are set to zero and the bosonic cou- 
plings are kept at the Standard Model values (fermio- 
phobic Higgs model). The largest excess with respect to 
the background-only hypothesis is found at 125.5 GeV, 
with a local significance of 2.9 standard deviations, 
which reduces to 1.6 standard deviations when taking 
into account the look-elsewhere effect. The data exclude 
the fermiophobic Higgs model in the ranges 110.0 - 
118.0 GeV and 119.5 - 121.0 GeV at 95% confidence 
level. 
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Several extensions of the Standard Model (SM) have 
been proposed in which the Higgs field couplings to 
some or all fermion generations are substantially sup- 
pressed, for example two Higgs doublet models or Higgs 
triplet models QUI] ■ A fermiophobic benchmark model, 
in which the Higgs field couplings to all fermions are 
set to zero while the couplings to bosons are kept at 
their SM values, has been introduced to allow a generic 
investigation of these scenarios [5]. 

In such a model, the production of the Higgs bo- 
son in hadron colliders and its decay properties are 
significantly altered compared to the SM. Fermiopho- 
bic Higgs bosons can only be produced via vector bo- 
son fusion (VBF) or associated production with vector 
bosons (VH, V = W,Z). Because Higgs boson decays 
to fermions are absent at tree level, the branching frac- 
tions for decays to gauge bosons are enhanced. In ad- 



dition, the partial width of the decay to two photons 
is enhanced by the suppression of the destructive in- 
terference between the VF-boson and top-quark loops. 
The resulting cross section times branching ratio for 
fermiophobic Higgs boson production with decay to two 
photons is larger than that of the SM for Higgs boson 
masses (run) below 125 GeV. Table [T] lists, for several 
values of uih , the fermiophobic Higgs boson cross sec- 
tion multiplied by the decay branching ratio into two 
photons. The ratio of this quantity with respect to that 
of the SM Higgs boson and the enhancement of the 
diphoton branching ratio are also shown. In addition to 
the enhanced diphoton decay rates, the recoiling jets or 
vector bosons in the VBF or VH production modes, re- 
spectively, imply a high transverse momentum for the 
Higgs boson that can be exploited as a discriminating 
variable in the analysis. However, for increasing m# the 
diphoton decay rate falls rapidly, making the search less 
sensitive at higher masses in this decay channel. 

Searches for a fermiophobic Higgs boson have been 
performed at the LEP and Tevatron colliders. The com- 
bination of results from the LEP experiments [5] ex- 
cludes a fermiophobic Higgs boson at 95% confidence 
level (CL) for masses below 109 GeV. When includ- 
ing both the WW and 77 decay modes, the Tevatron 
experiments exclude a fermiophobic Higgs boson with 
masses up to 119 GeV [6j[7]. 

This letter describes a search for a fermiophobic 
Higgs boson using diphoton events produced in proton- 
proton collisions at a centre-of-mass energy of y/s = 
7 TeV using data corresponding to an integrated lumi- 
nosity of 4.9 fb _1 collected by the ATLAS experiment. 
This analysis follows exactly that of the related search 
for a SM Higgs boson with the same dataset [5], but 
the fermiophobic Higgs hypothesis is used to construct 
the signal model. The sensitivity to the fermiophobic 
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Table 1: Higgs boson production cross section multi- 
plied by the branching ratio into two photons for the 
fermiophobic benchmark model (of), the ratio of this 
value to the SM value (<J{/<jsm) and the two pho- 
ton branching ratio enhancement compared to the SM 
(B{/Bsm) for various fermiophobic Higgs boson masses. 
The expected number of signal events after candidate 
selection are also shown for 4.9 fb -1 of data as well as 
the overall signal selection efficiencies. 



m H [GeV] 110 115 120 125 130 135 140 145 150 

cr f [fb] 163 90 53 32 21 13 8.9 5.9 3.9 

o-f/o-sM 3.7 2.1 1.2 0.8 0.6 0.4 0.3 0.3 0.2 

Bf/BsM 30.2 17.0 10.3 6.7 4.7 3.5 2.8 2.3 2.0 

Signal events 255 149 91 58 38 25 17 12 7.9 

Efficiency [%] 32 34 35 37 38 38 39 40 42 



signal is larger than that for the SM Higgs due to the 
larger diphoton transverse momentum. 

The ATLAS detector is described in detail in Ref. (9| . 
The most relevant subsystems for this analysis are the 
calorimeter, in particular the electromagnetic section, 
and the inner detector. The electromagnetic calorime- 
ter is a lead-liquid-argon detector, finely segmented 
in the lateral and longitudinal directions. It is com- 
posed of a barrel part covering the pseudorapidity range 
|?7| < 1.475 and two end-cap sections covering 1.375 < 
|?7| < 3.2. The barrel (I77I < 0.8) and extended barrel 
(0.8 < \rj\ < 1.7) hadron calorimeter sections consist of 
steel and scintillating tiles, while the end-cap sections 
(1.5 < |?/| < 3.2) are composed of copper and liquid 
argon. The inner detector includes silicon-based pixel 
and micro-strip detectors in the range | | < 2.5, and a 
transition radiation tracker with electron identification 
capability extending out to I77I < 2.0. It is surrounded 
by a superconducting solenoid that provides a 2 T axial 
magnetic field. 

Data used in this analysis were recorded using a 
diphoton trigger with a 20 GeV transverse energy (Et) 
threshold on each photon. This trigger is seeded by 
a first-level trigger, which requires two clusters in the 
electromagnetic calorimeter with Et > 14 GeV or Et > 
12 GeV, depending on the data-taking period. This trig- 
ger has a signal efficiency close to 99% following the 
final event selection. After application of data-quality 
requirements the analysed data sample corresponds to 
a total integrated luminosity of 4.9 ± 0.2 fb -1 [TU] , 

The events are required to have at least one recon- 
structed vertex with a minimum of three associated 



tracks, where the transverse momentum of each track 
is required to be larger than 0.4 GeV. At least two pho- 
tons within the fiducial region \rj\ < 2.37 (excluding 
the transition region between the barrel and the end- 
cap, 1.37 < \ri\ < 1.52) satisfying tight identification 
criteria based on electromagnetic shower shapes [TT] 
are required. The transverse momenta for the leading 
and sub-leading photons are required to be larger than 
40 GeV and 25 GeV, respectively. The photon recon- 
struction and identification efficiency ranges typically 
from 65% to 95% for Et in the range between 25 GeV 
and 80 GeV. The transverse energy deposited around 

each photon within a cone of AR = \J (Ar]) 2 + {Atfi) 2 
0.1. excluding the deposits of the photon itself, is re- 
quired to be less than 5 GeV. Corrections for the small 
estimated energy leakage outside the excluded region, 
the underlying event and effects of additional minimum 
bias interactions occuring in the same or neighbouring 
bunch crossings (in-time and out-of-time pileup) are ap- 
plied to this quantity on an event- by- event basis. 

The invariant mass of each diphoton candidate (/n 77 ) 
is evaluated using the photon energies, the impact points 
measured in the calorimeter and the production vertex. 
The photon energy calibration is performed indepen- 
dently for converted and unconverted photons. Con- 
verted photons are defined to be those with a well- 
reconstructed conversion vertex in the inner detector. 
A detailed simulation of the detector geometry and re- 
sponse is used for the calibration. Additional correc- 
tions due to mis-modelling of the material in front of 
the calorimeter and of calorimeter non-uniformities are 
applied. These amount to about ±1% depending on the 
pseudorapidity of the photon and are obtained from 
studies of Z — s- e+e~ decays in data [T2]. The dipho- 
ton production vertex along the beam axis is deter- 
mined by combining the trajectories of each photon, 
measured using the longitudinal segmentation of the 
calorimeter, with a constraint from the average beam 
spot position. The position of the conversion vertex is 
also used where the photons convert in the tracking 
region instrumented with silicon detectors. Conversion 
candidates with tracks reconstructed in inactive regions 
of the innermost pixel layer are rejected to reduce the 
contamination from misidentified electrons. The reso- 
lution of the diphoton mass reconstructed using this 
method is dominated by the photon energy resolution. 

A total of 22,489 events were selected with a dipho- 
ton invariant mass between 100 GeV and 160 GeV. Al- 
though not used directly in the final result, the diphoton 
sample composition was studied using a two-dimensional 
side-band technique based on photon identification qual- 
ity and isolation [5] . The fraction of true diphoton events 
was estimated to be (71 ± 5)%. The rest of the back- 



3 



ground is due to events with one or more misidenti- 
fied jets, except for a small (~ 0.7%) contribution from 
Drell-Yan events where both electrons pass the photon 
selection. 

To enhance the sensitivity of the analysis, the data 
sample is split into nine categories, each with differ- 
ent expected signal mass resolutions, signal yields and 
signal-to-background ratios (S/B). This categorisation 
depends on the impact point of the photons on the 
calorimeter, the presence of photon conversions and the 
value of the component of the diphoton transverse mo- 
mentum orthogonal to the diphoton thrust-like axis in 
the transverse planc0 (prt ) [SHU] ■ 

Events in which both photons are unconverted are 
separated into the unconverted central (both photons 
in the central region of the barrel calorimeter, |ry| < 
0.75) and unconverted rest (all other events) categories. 
Events for which at least one photon is converted are 
separated into the converted central (both photons within 
\i]\ < 0.75), converted transition (at least one photon 
close to the barrel/end-cap transition region, 1.3 < 
\r)\ < 1.75) and converted rest (the remaining events) 
categories. 

With the exception of the converted transition cat- 
egory, all the events are further subdivided into low pTt 
(PTt < 40 GeV) and high pTt (all other events) cate- 
gories. Monte Carlo (MC) simulation studies show that 
a fermiophobic Higgs boson signal has larger prt on av- 
erage than background events. This quantity is strongly 
correlated with the diphoton transverse momentum but 
offers several advantages. Higher values of pxt do not 
include kinematic configurations for which the two pho- 
tons are back-to-back in the azimuthal plane with sub- 
stantially different transverse momenta. This reduces 
biases on the identification and isolation (transverse en- 
ergy deposited around the photon) of the sub-leading 
photon in the high pxt categories and retains a mono- 
tonically falling diphoton invariant mass distribution 
for the background events at the chosen cut values. The 
latter quality is advantageous for the background mod- 
elling and associated uncertainties discussed below. 

A full Geant4-based [T3] MC simulation [I5J of Higgs 
boson events decaying into two photons is used to model 
the expected signal. The signal yields are normalised 
to next-to-next-to-leading-order production cross sec- 
tions [TTH22] and the branching ratios for the fermio- 
phobic Higgs boson are calculated using HDECAY |23J. 
Higgs boson VBF production is simulated using POWHEG [24J 
interfaced with PYTHIA [25] for showering and hadro- 

1 PTt = lp?> 7 x tli where t = — ^rrr denotes the transverse 

11 1 Ipt — Pt I 

thrust, p T ' and p T 2 are the transverse momenta of the two 
photons, and p T 7 = p T ' + p T 2 is the transverse momentum 
of the diphoton system. 



nisation, while PYTHIA is chosen for the VH processes. 
Pileup effects are simulated by overlaying each MC event 
with a variable number of simulated inelastic pp colli- 
sions, taking into account the LHC bunch-train struc- 
ture [25]. 

A set of corrections is applied to the simulated events 
in order to match the data-taking conditions. The sim- 
ulated events are re-weighted to reproduce the distri- 
bution of the average number of interactions per bunch 
crossing reconstructed in the data, which has a mean 
value of about nine for the data sample used in this 
analysis. The energies of the simulated photons are smeared 
to account for differences observed in studies of the 
calorimeter resolution with Z — > e + e~ decays. Calorime- 
ter shower shapes used in the photon identification are 
slightly shifted to improve the agreement with the dis- 
tributions observed with inclusive photons from data. 

The number of fermiophobic Higgs bosons expected 
after candidate selection and the overall signal selection 
efficiency for various values of mu are shown in Table [T] 
The signal selection efficiency increases from 32% to 
42% as the Higgs boson mass increases from 110 GeV 
to 150 GeV. 

Table 2: Expected signal mass resolution (ctcb and 
FWHM in GeV, see text) and total number of sig- 
nal events (iVg) for win = 120 GeV for each of the 
nine analysis categories and for the inclusive case. Also 
shown for each category are the number of observed 
events (N-a) m the diphoton mass range from 100 GeV 
to 160 GeV, and the expected signal-to-background ra- 
tio (S/B) in a mass window containing 90% of the sig- 
nal. 



Category 


"CB 


FWHM 


N s 


N u 


S/B 


Unconverted central, low pxt 


1.4 


3.3 


6.2 


1763 


0.03 


Unconverted central, high pxt 1-3 


3.2 


8.6 


235 


0.37 


Unconverted rest, low pxt 


1.7 


3.9 


12.1 


6234 


0.02 


Unconverted rest, high pxt 


1.6 


3.8 


16.0 


1006 


0.13 


Converted central, low pxt 


1.6 


3.8 


4.0 


1318 


0.02 


Converted central, high pxt 


1.5 


3.5 


5.8 


184 


0.26 


Converted rest, low pxt 


2.0 


4.6 


11.8 


7311 


0.01 


Converted rest, high pxt 


1.9 


4.4 


16.1 


1072 


0.09 


Converted transition 


2.3 


5.8 


10.8 


3366 


0.01 


All categories 


1.7 


3.9 


91.2 


22489 


0.03 



The signal is modelled as the sum of a core compo- 
nent, described by a Crystal Ball (CB) function [27], 
and a wider Gaussian component incorporating outly- 
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ing events. The latter component typically accounts for 
less than 5% of the signal. Table[2]lists the expected full- 
width-at-half-maximum (FWHM) and Gaussian width 
of the core component (ctcb) for each of the nine event 
categories. The expected number of signal events for 
ma = 120 GcV, the number of background events in 
the diphoton mass range of 100 GeV to 160 GeV, and 
the signal-to-background ratio in a mass window con- 
taining 90% of the signal are also shown. The main sen- 
sitivity to the fermiophobic production modes comes 
from the high prt categories due to their enhanced 
signal yields and signal-to-background ratios. Figure Q] 
shows the signal diphoton mass distribution summed 
over the high pxt categories for a Higgs boson mass of 
120 GeV. 
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Fig. 1: Diphoton invariant mass spectrum from sim- 
ulated signal samples (dots) with m# = 120 GeV 
summed over the high pit categories, superimposed 
with the signal model (line). 



The observed diphoton invariant mass distribution 
in each category is modelled by an exponential func- 
tion. A fit to the data is performed for which the slope 
and normalisation are unconstrained. Studies with large 
samples of simulated diphoton events show that this 
simple function gives a good description of the expected 
shape. The small systematic uncertainties associated 
with this assumption are discussed below. Figures [U[a) 
and Eljb) show the diphoton mass distributions of the 
selected data events summed over the low and high pxt 
categories, respectively. The converted transition cate- 
gory is included in the low prt categories. 

Systematic uncertainties affecting the signal signif- 
icance arise from uncertainties on the predicted signal 
yields, the expected partition of the signal among the 
categories and the modelling of the signal and back- 
ground shapes. The dominant experimental uncertainty 
on the signal yield is due to the imperfect knowledge 
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Fig. 2: Diphoton invariant mass spectra for the (a) low 
and (b) high pxt categories, overlaid with the sum of 
the background-only fits from the individual categories. 
The bottom plots show the residual of the data with 
respect to the fitted background. The signal expectation 
for a fermiophobic Higgs boson with a mass of 120 GeV 
superimposed on the background fit is also shown. 



of the photon reconstruction and identification efficien- 
cies, which is estimated to be ±11%. This uncertainty 
is studied with electrons from W and Z boson decays 
in data, and photons from radiative decays of Z bosons 
to electron and muon pairs. In addition, the effect of 
pileup on photon identification gives a further contri- 
bution to the signal yield uncertainty of ±4%. Uncer- 
tainties related to the trigger efficiency (±1%), isolation 
cut efficiency (±5%) and luminosity (±3.9%) are also 
included here. 

Uncertainties on the signal cross section include a 
combination of the uncertainties on the parton distribu- 
tion functions 28, 25] and a s , and uncertainties on the 
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QCD scale. Combining the VBF and VH production 
modes this uncertainty is within ±4% over the con- 
sidered mass range. To this uncertainty, that due to 
the H — > 77 branching ratio (±5%) is added linearly, 
based on the SM calculation 22 . This yields uncertain- 
ties of ±9% on the theoretical signal yield, leading to an 
overall uncertainty of ±16% on the total signal expec- 
tation. In addition, the uncertainty on the Higgs boson 
Pt modelling is estimated by comparing signal samples 
from alternative MC generators - HERWIG [30] for VBF 
and ResBos [3T] for VH. The result is a ±1% signal mi- 
gration between the low and high pTt categories with a 
negligible effect on the signal selection efficiency. 

The dominant uncertainties on the signal mass res- 
olution are due to the uncertainty on the calorime- 
ter energy resolution (±12%) and photon calibration 
(±6%), which are both extrapolated from the uncer- 
tainty on the electron calibration determined using Z 
and J/4> data [12J. The latter comes from the imper- 
fect knowledge of the material in front of the active 
part of the calorimeter and is estimated using simula- 
tions with different amounts of material. This quantity 
also affects the fraction of expected events in the cate- 
gories with converted photons; the maximal migration 
between converted and unconverted categories is esti- 
mated to be ±4.5%. Other effects on the signal mass 
resolution are due to pileup fluctuations contributing 
to the cluster energy measurement (±3%) and to the 
uncertainty on the photon angular resolution (±1%) 
which is studied in Z — > e + e~ decays using the track- 
based direction measurement. The total relative uncer- 
tainty on the diphoton invariant mass resolution is thus 
±14%. 

Systematic uncertainties on the background mod- 
elling arise from a possible deviation of the background 
mass distribution from the assumed exponential shape. 
This uncertainty is evaluated as the number of events 
that could be mistakenly attributed to the signal. It is 
estimated from the adequacy of the chosen background 
model's description of the mass distribution predicted 
by ResBos [35]. The residuals of the fit of the back- 
ground model to the ResBos diphoton mass distribution 
are integrated over a sliding mass window of 4 GeV, the 
approximate FWHM of the expected signal. The largest 
deviations were found at small invariant masses and 
these uncertainties are then applied over the whole mass 
range. The resulting uncertainties range from ±0.1 to 
±7.9 events in the individual analysis categories, where 
the magnitude of these uncertainties is roughly propor- 
tional to the number of background events in each cate- 
gory. These absolute uncertainties do not scale with the 
signal strength in the final likelihood fit. For a fermio- 
phobic Higgs boson with win — 120 GeV the back- 



ground modelling uncertainty in the high pTt categories 
is equivalent to up to 5% of the signal yield with nomi- 
nal signal strength. The estimation of the uncertainties 
is cross-checked by fitting the data with different func- 
tional forms and comparing the result to the exponen- 
tial fit. 

The possible presence of a signal is investigated us- 
ing a combined likelihood function constructed from the 
signal and background models for the diphoton invari- 
ant mass distribution in each of the nine categories. Un- 
binncd maximum likelihood fits of the signal strength 
are performed, treating the systematic uncertainties as 
nuisance parameters - fourteen in total. These nuisance 
parameters are added to the signal likelihood function 
using a Gaussian term for the background modelling 
uncertainty, and log-normal terms for all other uncer- 
tainties. 
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10 f ; ■ Obs. p (Low p T( + Conv. transition) ^ 
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m H [GeV] 

Fig. 3: Local observed po as a function of the Higgs 
boson mass mn (solid line) and the median expectation 
for a fermiophobic signal with the given run (dotted 
line). The five points near 125 GeV show the observed 
Po when the uncertainty on the photon energy scale 
is considered. The individual contributions of the low 
PTt and high pTt categories to the observed po are also 
shown. 



The compatibility of the data with the background- 
only hypothesis, relative to the hypothesis of background 
plus the fermiophobic model signal, is quantified by the 
local significance po- Figure [3] shows the result for mu 
ranging from 110 GeV to 150 GeV, where po is com- 
puted in 0.5 GeV steps using asymptotic formulae [3*3] . 
The contributions to po values from the high pTt and 
low pxt categories are shown separately. The high pTt 
contribution has a minimum po at 125 GeV, while the 
low pxt contribution has a minimum at 127 GeV. The 
larger signal-to-background ratio as well as the larger 
expected signal yield in the high pTt category compared 
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to the low pTt category results in the high pxt contribu- 
tion dominating in the final result. The combined po has 
a minimum at 125.5 GeV corresponding to 3.0 standard 
deviations. The figure also shows the po value expected 
for a fermiophobic Higgs boson signal, as a function of 
Higgs boson mass. 

To obtain the final result, the impact of the un- 
certainties on the photon energy scale is considered for 
Higgs boson masses in the region of the minimum po, as 
shown in Fig. [3] The corresponding effect on the mea- 
sured po value is estimated using pseudo-experiments, 
since asymptotic formulae were found not to yield ac- 
curate estimates of the probability in this case. The 
position of the minimum po is almost unchanged and 
the significance is lowered to 2.9 standard deviations. 
Taking the look-elsewhere effect [33] into account in 
the range 110 - 150 GeV, the significance reduces to 
about 1.6 standard deviations, with po ~ 0.051. This 
may be compared to the result of a search for the SM 
Higgs boson performed with the same dataset and can- 
didate selection [5], yielding a minimum pq at a mass 
of 126.5 GeV with a global significance of 1.5 standard 
deviations. No statistically significant preference for ei- 
ther the SM or fermiophobic Higgs boson is observed. 

Given the lack of evidence for a signal, mass-dependent 
exclusion limits on the fermiophobic benchmark model 
are calculated at the 95% confidence level (CL) with a 
profile likelihood ratio test statistic in the CL S modified 
frequentist approach 33,35, 36 and are shown in Fig. 2J 
Fermiophobic Higgs boson masses from 110.0 GeV to 
118.0 GeV and from 119.5 GeV to 121.0 GeV are ex- 
cluded, while the expected exclusion mass range is 110.0 
- 123.5 GeV. These results give more stringent lower 
mass limits than the previous results from LEP 
(108.2 GeV) [5] and the Tevatron (112.9 GeV from DO, 
114 GeV from CDF) [61137] in the diphoton decay chan- 
nel. 
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